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ABSTRACT 

In the summer of 2012, during a Pulsar Search Collaboratory workshop, two high-school students discovered 
J1930-1852, a pulsar in a double neutron star (DNS) system. Most DNS systems are characterized by short 
orbital periods, rapid spin periods and eccentric orbits. However, J1930-1852 has the longest spin period 
(Apin ~185 ms) and orbital period {P^ ^45 days) yet measured among known, recycled pulsars in DNS systems, 
implying a shorter than average and/or inefficient recycling period before its companion went supernova. We 
measure the relativistic advance of periastron for J1930-1852, a) = 0.00078(4) deg/yr, which implies a total 
mass (Mtot = 2.59(4) M©) consistent with other DNS systems. The 2a constraints on Mtot place limits on the 
pulsar and companion masses (iWp < 1.32M0 and me > 1.3OM0 respectively). J1930-1852’s spin and orbital 
parameters challenge current DNS population models and make J1930-1852 an important system for further 
investigation. 


1. INTRODUCTION 

To date, ^2,300 pulsars are known (|Hobbs et al.|2004jl and 
^10% of them are in binary systems, orbiting white dwarf 
(WD), neutron star (NS) or main sequence star (MS) com¬ 
panions. The vast majority of these binaries are NS-WD sys¬ 
tems; many of these systems emerge from scenarios where 
the pulsar forms first, followed by its companion, which over¬ 
flows its Roche Lobe; accretion transfers angular momen¬ 
tum to the pulsar, decreasing the spin period a nd resulting 
in a millisecond pulsar (MSP) orbiting a WD (|Alpar et al. 
1982|l. This process of accretion and spin-up is commonly 


referred to as recycling and the period derivative of a recy¬ 
cled pulsar tends to be significantly lower than that of an un¬ 
recycled pulsar with the same spin period. There are four 


evolve off the main sequence Johnston et al. 

1992 |Kaspi| 

et al.|1994l Stairs et al.|2001| Lyne|20051l; an ad 

ditional four 

have been found with planet-sized companions (|1 

’horsett et al.| 

|19931|Wolszczan|1994||Bailes et al.|201 IHStoval 

1 et al.|2014|l. 

More massive companions end their evolution of 

' the main se- 


quence in supernovae, resulting in double neutron star (DNS) 
systems. DNS systems are far less likely to remain bound than 
NS-WD systems, since the former must survive two super¬ 
nova explosions during formation. Only about 10% of these 
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binary system s remain bound after one supernova explosion 
( Bailes|1989 1. The probability of remaining bound after two 
supernovae is much lower (^1%) and only nine such systems 
have been found and studied previously (see references in Ta¬ 
ble 

DNS systems have tantalizing applications — for exam¬ 
ple, testing theories of gravity by measuring relativistic ef¬ 


fects ( [Fonseca et al.|20141 Weisberg et al.|2010|[Kramer et al. 

and predicting DNS merger rates relevant to ground- 

gravitational wave detectors like LIGO (|Kim et ak] 


2010| 2013[). DNS systems have also provided some of the 


most precise NS mass measurements, allowing for a statistical 
investigation of the underlying mass distribution ( Ozel et al] 


2012[ Sdiwab et al.|2010[|'lTiorsett & Chakrabartyjl99^i. 


Thorsett & Chakrabarty| ( T^99| l used a sample of 26 NSs (21 

MSRs and hve binary companion NSs) with measured masses 
to deter mine a mean NS mass , (m) = L35±O.O4M0. More 
recently, [Schwab et al.| (|2010]l argue that the underlying NS 
mass distribution for objects in DNS systems is bimodal, with 
narrow peaks at 1.246 M0 and 1.345 M0. They also suggest 
that these peaks indicate unique formation scenarios, where 
the lower mass c omponent represents NSs that formed via 
electron capture ( |Nomoto||1984[ |Podsiadlowski et ^|2004| l 
and the higher mass component is i ndicative of iron c ore- 


Ozel et al. 


(20121 use 


collapse ( [Woosley & Weaver) 1986[ ). ^^_ 

a Bayesian statistical approach to infer mass distributions for 
NSs with distinct evolutionary histories; they agree that NS 
masses provide clues about respective formation scenarios, 
however, they e xpress skepticism that NS mass distributions 
Schwab et al.| ( |2OT0] l claim. Therefore, ad- 


are as narrow as 


ditional information is necessary in identifying a NS’s evolu- 
tionary history. 


I Wong et al. (|2010]l investigate core-collapse mechanisms 
in eight Galactic DNS systems by inferring progenitor mass 
of the second-bom NS and the magnitude of the supernova 
kick it received at birth from measured DNS orbital pa¬ 
rameters and kinematic information. Using these methods, 
they conclude that NS companions of PSRs B1534H-12 and 
B19134-16 underwent iron core-collapse supernovae, while 
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J0737-3039A’s companion likely formed via electro n cap- 
ture supernova. This final result was corroborated by |Ferd- 
man et al. (2013[); through detailed pulse profile shape analy¬ 


sis, they constrained the double pulsar system geometry, con¬ 
cluding that the secondary supernova explosion was relatively 
symmetric, indicative of an electron capture process. 

There is a long history of work contributing to the ide a that 
J0737-3039B formed via electron capture supernova. |Pod- 
T] ( 2004| l first suggested that the critical stel- 
d to foiT 


siadlowski et al 
lar mass requirei 


Toim a NS (10-12 Mq for solitary stars) 
should be significantly lower for tight, interacting binary sys¬ 
tems (6-8 M 0 ). Based on early scintillation velocity mea¬ 
surements, contraints were placed on the progenito r mass of 
J0737-3039B and kick velocity due to its supernova ([Willems 


& Kalogera|2004| Piran & Shaviv|20()5 Willems et al.|200^ 


2006|. Precise transverse velocity measurements from an ex 

tended timing campaign provided the necessary information 
to claim an unusually low progenitor mass for J 0737-3039B 
and corresponding low supernova kick velocity (jStairs et al. 
[2006l|Piran & Shaviv|2006) . 

We draw attention to the double pulsar system here to il¬ 
lustrate the detailed process that is required to make claims 
about DNS formation scenarios. As a result of the mass con¬ 
straints presented in this paper, 11930-1852 appears to be in a 
DNS system. Its unique spin and orbital parameters challenge 
models that describe DNS formation. 

In ® we describe the GBT 350 MHz Drift Scan survey and 
the Pulsar Search Collaboratory, as well as the follow-up tim¬ 
ing campaign and the parameters measured for J1930-1852; 
Oprovides evidence that a NS companion is likely, although 
rSio follow-up observations have not provided any evidence 
of a pulsar companion. In ^ we draw conclusions from our 
findings and outline plans for future work. 

2. TIMING OBSERVATIONS & ANALYSIS 

In May-August of 2007, when the Green Bank Telescope 
(GBT) was undergoing trac k replacement, the GBT 350-MHz 
Drift Scan Pulsar Survey (|Boyles et al.||2013l [Lynch et al. 
2013|l looked for radio pulses as the sky drifted overhead. Of 
the 1,491 hours of recorded drift scan data.~300 hours were 
allocated to the Pulsar Search Collaborator}j^(PSC). The sur¬ 
vey and follow-up timing observing campaigns, processing 
pipeline an d the students’ first five discoveries are discussed 
in detail in |Rosen et al. 1 ( 2013) 1. PSR J1930-1852 is the sixth 
pulsar discovered by PSC students. 

After J1930-1852 was flagged by PSC students in late July 
of 2012, we confirmed the candidate with a long scan at 
350 MHz with the GBT, at which point we also verified the 
binary nature of the source, given the significantly different 
measured spin period than that reported on the discovery plot. 
Soon after confirmation, we performed “gridding" observa¬ 
tions, tiling the 350 MHz beam with seven 820 MHz beam po¬ 
sitions, since the higher-frequency receiver has a smaller an¬ 
gular beam size on the sky (^0.25 de g compared to ~0.5 deg). 
As described by|Morris et al.|(|2002|l, gridding reduces the un¬ 
certainty on the pulsar’s position and makes it easier to even¬ 
tually achieve a phase-connected timing solution. At the GBT, 
the 350 MHz prime-focus receiver is usually only mounted for 
several days each month, while the 820 MHz receiver is up for 



Fig. 1.— Black points show spin period measurements at various observing 
epochs with error bars, often smaller than the points themselves. The gray 
dashed line illustrates predicted spin period versus time based on our five 
Keplerian, orbital parameters. The apparent spin period variation shown here 
is a Doppler effect due to binary motion and the pulsar’s intrinsic spin period 
is given in Table^ 


TABLE 1 

Timing parameters for PSR J 1930-1852 


Measured Parameters 

Value 

Right Ascension (J2000). 

Declination (J2000). 

Spin Period (s). 

Period Derivative (s/s). 

Dispersion Measure (pc cm“^). 

Reference Epoch (MJD). 

Span of Timing Data (MJD). 

Number of TOAs. 

RMS Residual (fis) . 

xL . 

19:30:29.7156(7) 

-18:51:46.27(6) 

0.18552016047926(8) 

1.8001 (6) X 10^'^ 
42.8526(4) 

56513 

56121-56904 

75 

29 

1.05 

Binary Parameters 


Orbital Period (days). 

Projected Semi-major Axis (It-s) ... 

Epoch of Periastron (MJD). 

Longitude of Periastron (deg). 

Orbital Eccentricity. 

Advance of Periastron (deg/yr). 

45.0600007(5) 

86.890277(7) 

56526.642330(3) 

292.07706(2) 

0.39886340(17) 

0.00078(4) 

Derived Parameters 


Surface Magnetic Field (10^® Gauss) 
Spin-down Luminosity (10^^ erg/s). 

Characteristic Age (Myr). 

Mass Function (Mq). 

Minimum Companion Mass (Mq) .. 

Combined Mass (Mq) . 

Mean S 820 (mjy). 

6.0 

1.1 

163 

0.34690765(8) 

IS, 

0.7 


Note. — Uncertainties in the last significant digit(s) are quoted in 
parentheses and represent Icr en'ors on measured parameters. Since 
the flux density (5820) quoted here is based on a single-epoch mea¬ 
surement, the uncertainty may be 10-20%. 

^Minimum comp anion mass listed here is based on constraints pro¬ 
vided by Mtot (see ^3.2[ l. 


The PSC jRosen etal .|2010| aims to interest high-school students in sci¬ 
ence, technology, engineering and mathematics (STEM) related career paths, 
focusing especially on engaging women and minority students as well as 
those from low-income families. 


the rest of the month, so gridding also provided more flexibil¬ 
ity in our follow-up timing campaign. 

We first conducted high-cadence, then monthly timing ob¬ 
servations once we had an orbital solution. We observed 
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M.jD-56513 (days) 

Fig. 2.— Timing residuals for J1930—1852 plotted here correspond to fit parameters listed in Table[T] Red and black points represent 350MFlz and 820MHz 
observations respectively. The group of 350 MHz TOAs at MJD 56904 come from a 3.5 hr observation at superior conjunction. 


J1930-1852 at 350 MHz and (primarily) 820 MHz center fre¬ 
quencies, with 100 MHz and 200 MHz of bandwidth respec¬ 
tively. For all observations, we used th e Green Bank Ultimat e 
Pulsar Processing Instrument (GUPPI; |DuPlain et al.|(|2008|l) 
with 81.92 fis resolution time and 2048 frequency channels. 

During each session, we observed J1930-1852 for 
~ 15 mins, manually excised RFI with psrzap (part o f the 
PSRCHIVlj^ software package, Hotan et al. (2004l) and 


then summed the signal across the entire bandwidth. We 
summed across the time domain to generate one mean pulse 
profile per session and compared it with a synthetic stan¬ 
dard profile to compute a time of arrival (TOA) using the 
PSRCHIVE routine pat. Standard profiles, one for each 
observing frequency, were created by htting Gaussian com¬ 
ponents to a high signal-to-noise prohle. 

We measured a spi n period at each observing epo c h then , 
using methods from |Bhattacharyya & Nityananda] (|2008|l, 
found an orbital solution for J1930-1852; period rneasure- 
ments and those predicted by our orbital solution are shown in 
Figure [T] We furthe r rehned orbital parameters in TEMP 02 p] 
( |Hobbs et al. |2006| l to achieve a full, phase-connected solu- 
tion, shown in Table H Our timing solution, in B ary cen¬ 
tric Coordinate Time ^CB), uses the DE405 Solar System 
ephemeris and TT(BIPM) clock corrections. Subtracting 
TOAs from modeled arrival times determined by our phase- 
connected timing solution yields timing residuals shown in 
Figure]^ Since the reduced chi-squared statistic, Xred ~ ^ 
our timing residuals, we do not use a multiplicative “error fac¬ 
tor” (EFAC), so uncertainties given in Table[T]are identical to 
those reported by Tempo2. 

On September 4, 2014, we observed J1930-1852 at supe¬ 
rior conjunction (orbital phase, (f) ~ 0.373) for 3.5 hrs with our 
normal 350 MHz setup described above. We obtained several 
TOAs from this epoch - displayed in Eigure and included 
in the timing solution reported in Table Since the post-fit 
residuals did not show any sign of a Shapiro delay signature, 
we do not fit for Shapiro delay range (r) and shape {s) param¬ 
eters here. 

A flux- and polarization-calibrated, 820-MHz pulse pro¬ 
file for J1930-1852 is shown in Eigure About 15% of 


http://psrchive.sourceforge.net 

http://tempo2.sourceforge.net 



Pig. 3.— Bottom panel: A coherently dedispersed, flux- and polarization- 
calibrated pulse profile for J1930—1852 obtained from a 1 hr observation at 
820 MHz with 200 MHz bandwidth and 1024 profile bins. The black line 
represents total intensity, while linear and circular polarization are shown in 
blue and red respectively. Top panel: The position angle (P.A.) swing due to 
slight linear polarization in the leading edge of the pulse profile. 


J1930-1852’s emission is linearly polarized, but due to the 
relatively flat position angle curve, we w ere not able to fit 
this curve using the rotating vector model (|Radhakrishnan & 
Cooke|1969| l. There is no trace of circularly-polarized emis- 
sion. 


3. NATURE OP THE COMPANION 

Given optical images of the sky surrounding J1930-1852 
and mass constraints based on our timing solution (Table [T]), 
the pulsar’s companion is most likely another neutron star. 

3.1. Optical Follow-up 

Assuming the companion is a main sequence star, we 
can estimate its apparent bolometric magnitude and appar¬ 
ent SDSS g magnitude. Eirst, with the mass-luminosity rela- 
tion for an ap propriate mass range Lc.min/L© = (nic,min/M 0 )^-^ 
(|Allen|l973|l, where mc.min is the minimum companion mass, 
we hnd Lc.min = 2.5 L©, which would correspond to an E5V 
spectral type. Next, using a DM-estimate distance, (7 dm ~ 
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Fig. 4.— This 1.9'X 2.5' field of view sun'ounding J1930-1852 is the result 
of stacking 20 300 s images taken with a SDSS g filter and the CTIO 0.9 m 
telescope. No objects were detected within a five-arcsecond radius around 
J1930—1852 — well beyond the uncertainties in measured position. 


1.5 kpc ( Cordes & Lazio[ 
metric magnitude of mboi., 


we estimate an apparent bolo- 
4.7. To convert this to an ap¬ 


parent g magnitude, we (i) applied a bolometric correction 
(BC ^ -0.09) to co nvert mtoi.c to an apparent V magnitude 
(Bessell et al. 1998 1, (ii) transformed this V magnitude to a 
SPSS g magnitude ( Jester et al.|2005 1, and (iii) included the 
effects of extinction (Ag ^ 0.33 ma^ using the estimates of 


Schlafly & Finkbeiner (20111. If the companion is a main se- 
quence star, we hnd that it should have an apparent magnitude 
of gc = 15.27 or brighter. 

We observed the held around J1930-1852 on the night of 
9 May 2014 using the CTIO 0.9 m telescope, w hich w as ac¬ 
cessed through the SMARTS ( jSuhasavage et al.||2010[) Con¬ 
sortium. We obtained 20 images of the held through a SPSS g 
hlter (“CTIO 4770/1006”) over the course of two hours; each 
individual exposure had an exposure time of 300 seconds, giv¬ 
ing us a total integration time of 6000 seconds. To minimize 
processing time, we read out only a 291 x 375 pixel subsection 
of the full 2048 X 2048 pixel CCP, which, with a 0.401 "/pixel 
plate scale, gave us a 1.9'x2.5' held of view (see Figure®. 

All frames were hat-h elded and bi as-subtracted using stan¬ 
dard routines in IRAlf®( Tody|l986 1, and the 20 reduced ob¬ 


ject frames were averaged together into a master frame with 
the imcombine task. We used astrometry.net to perform an 
astrometric calibration of the master frame to a precision of 
better than 0.1". To place an upper limit on the magnitude 
of any optical signatures, we hrst determined the aperture di¬ 
ameter that maximized the signal-to-noise (S/N) ratio in the 
photometry of nearby stars (2.4", or about two times the see¬ 
ing); next, we calculated the number of counts in the area of 


** IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As¬ 
tronomy, Inc., under cooperative agreement with the National Science Foun¬ 
dation. 


sky centered on 11930-1852 and determined the level of noise 
at this position. Setting a 5a detection threshold requirement 
for faint stars in the vicinity, we hnd a limiting magnitude of 
g = 20.5 for a possible optical counterpart of 11930-1852. 

Since the estimated g magnitude for a MS star companion 
gc = 15.27 is ~ 125 times brighter than the limiting magnitude 
achieved here, we rule out a MS companion. The uncertainty 
in doM 25%) does not alter this conclusion. 


3.2. Mass Constraints 

The mass function expresses the mass of the pulsar (Wp) 
and that of the companion (m^) in terms of Keplerian orbital 
parameters UpSin/ (projected semi-major axis) and P\, (orbital 
period). Using those measured parameters and setting the in¬ 
clination angle i = 90°, we place a lower limit on for any 
given iWp. 

The measurement of the relativistic advance of periastron 
(w) also provides the system’s total mass, = 2.59(4) Mq. 
Taking i = 90°, we use the mass function and double the un¬ 
certainties on Mtot to place 2cr lower/upper limits on the mass 
of the companion (me > I.3OM0) and the mass of the pulsar 
(mp < I.32M0) respectively. 

Based on the companion mass lower limit iWc ll = 1 -30 M 0, 
the orbital eccentricity e ~ 0.4, the spin period and period 
derivative of J1930-1852, a NS is the most likely companion. 


3.3. Radio Follow-up 

We dedispersed our 3.5 hr superior conjunction obser¬ 
vation using the DM measured for J1930-1852 (DM = 
42.85 pc cm“^), took a discrete Fourier Transform of the 
resulting timeseries and performed an acceleration search 
for a possible pulsar companion. To do so, we used the 
accelsearch routine from the PrestcF] suite of pulsar 
search software ( [Ransom 200 Ij ), examining me frequency do¬ 
main with a matched hlter template up to 20 Fourier bins 
wide. To maximize our sensitivity to pulsars with narrow pro- 
hles, we summed up to 16 harmonics for candidate signals. 
This procedure returned over 800 high-signihcance (S/N > 9) 
candidates, from which we removed obvious RFI and those 
that were harmonically related to J1930-1852 or each other. 

We folded and visually inspected the remaining candidates, 
but did not hnd any evidence of a pulsar counterpart. Assum¬ 
ing the harmonic summing was close to ideal, we were sensi¬ 
tive to S/N > 9 signals in the time domain, which corresponds 
to a 350 MHz Hux limit of ^30 pJy, given our observing set¬ 
up. These results suggest that a possible pulsar companion 
is either too weak to be detected or is not beaming along our 
line of sight. Given the sensitivity limit reached, the latter 
explanation is more likely. 


4. SUMMARY & CONCLUSIONS 

Analysis presented in ^implies that J1930-1852’s most 
likely counterpart is another NS. We compare J1930-1852’s 
parameters with those of other pulsars in known DNS 
systems that have similar evolutionary scenarios in Table 
1^ Given J1930-1852’s moderately short spin period and 
period derivative well below those measured for otherwise 
similar unrecycled pulsars, J1930-1852 is partially recycled 
and therefore, most likely formed before its companion. 
On a P — P diagram (see Figure®, J1930-1852 falls in the 
same region as other recycled DNS systems. Meanwhile, 


https://github.com/scottransom/presto 
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TABLE 2 

Known Recycled Pulsars in DNS Systems 


Pulsar 

^spin 

(ms) 

P 

(10-*** s/s) 

e 

^orb 

(days) 

Mp 

(M©) 

nic 

(Mq) 

A^tot 

(Mq) 


Recent 

References 



J0737-3039A 

22.7 

1.8 

0.09 

0.10 

1.34 

1.25 

2.59 


Kramer et al.l 

2006 

■ 


J1756-2251 

28.5 

1.0 

0.18 

0.32 

1.34 

1.23 

2.57 


Perdman et al. 

2014 

i 


B1913+16 

59.0 

8.6 

0.62 

0.32 

1.44 

1.39 

2.83 


weisbers et al. 



B1534-H12 

37.9 

2.4 

0.27 

0.42 

1.33 

1.35 

2.68 


Fonseca et al. 

zDTTf 

L-. 


J1829+2456 

41.0 

0.05 

0.14 

1.18 

< 1.34 

> 1.26 

2.53 

L.nampion et ai. |i 

01Wf005| 


J0453+1559 

45.8 

0.19 

0.11 

4.07 

1.54 

1.19 

2.73 

Martinez et al. (m pre 

T1 


J1518+4904 

40.9 

0.03 

0.25 

8.63 

< 1.17 

> 1.55 

2.72 


iJanssen et al.l 

2008 



J1753-2240 

95.1 

1.0 

0.30 

13.6 

- 

- 

- 


iKeith et al.l ^ 




J1811-1736 

104 

0.9 

0.83 

18.8 

< 1.64 

>0.93 

2.57 

1 

Lorongiu et al. 


\ 


J1930-1852 

186 

18 

0.40 

45.1 

< 1.32 

> 1.30 

2.59 


— 




Note. — A comparison between J1930—1852 and all other known primary, partially-recycled DNS pulsars, sorted by Pm-h. PSRs J1906+ 0746 and 
J0737—3039B were omitted because neither underwent recycling ^Lorimer et aL|20061|Kramer et al.|2006). PSRB2127+11C jJacoby et al.|2006) was also 
omitted because it was formed in a globular cluster, indicating a dilterent evolutionary history. 
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Fig. 5. — A P-P diagram s howing all pulsars in DNS systems (stars/squares) and all other known pulsars (dots). Measured P and P come from the ATNF Pulsar 
Catalog jHohbs et al.|2004f and lines of characteristic age and surface magnetic field are shown with dot-dash and dashed lines, respectively. Recycled DNS 
pulsars (stars) appear between the normal and millisecond pulsar populations and are listed in Table|^ Despite its significantly longer spin period, J1930-1852 
clearly belongs in the population of recycled DNS pulsars, unlike J1906-H0746 and J0737—3039B (squares) - neither of which have undergone recycling. 


it has a longer spin period (Pspin "^185 ms) and higher rate 
of spin-down (P ~ 2 x 10“'^ s/s) than any other first-born, 
recycled DNS pulsar and an orbital period (Pb ^ 45 days) 
longer than any other DNS system. These together may imply 
a shorter than average and/or inefficient mass transfer phase 
before the companion went supernova. Since accretion is 
thou ght to be the source of significant mass gain for recycled 
NSs ( Zhang et aLpOll) , a short and/or inefficient accretion 
period is consistent with the relatively low upper limit we 
place on the mass of J1930-1852, < 1.32M0. The upper 


mass limit for J1930-1852 is lower than those of all other 
recycled DNS pulsars, except J1518-1-4904. 


4.1. Spin Period-Eccentricity Relationship 


The Pspin - e r elationship was first n oted in McLaughlin 
let al.| (|2005() and |Faulkner et al.| (|2005|), after the discover- 
leTof J1756-2251 and J18293-2456; |Faulkner et dl ( |2005] l 
performed a linear fit in spin period/eccentricity space (see 
Figure!^, using the first seven known primary, partially recy¬ 
cled DNS pulsars and found a Pearson correlation coefficient 
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exhibits a slope simil ar to that of the empirical fipi n ~ g rela- 
tionship first noted by [McLaughlin et ah ( |2005| ) and Faulkner 
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Fig. 6. — The dashed line shows the spin period-eccentrici ty relation first 
noted by|McLaughlin et al. 1(2005) andjFaulkner et al. 1120051, which is sim¬ 
ply a fit using the first seven known primary, partially-recyclM DNS pulsars 
(unmarked points). The shaded region shows t he approximate sp read in prop¬ 
erties of simulated DNS systems according to|Dewi et al.|(20051, assuming a 
small Maxwellian kick velocity dispersion induced by the second supernova 
( 0-2 = 20km s-‘)0 PSRs J1753-2240 and J1930-1852 have been labeled 
because they bothiall off the dashed line, but have spin periods a nd eccen¬ 
tricities that are still roughly consistent with the|Dewi et al.|(20051 synthetic 
population. PSR J0453-H1559, a recently discovered DNS with a large mass 
asymmetry, described by Martinez et al. (in prep.), is consistent with the 
shaded region as well as the distribution of previously-discovered DNS pul¬ 
sars shown here. 

“Although we were unable to obtain actual population synthesis results from 
Figure 3 in that paper, we estimated the spread in spin period and eccentricity 
for simulated DNS systems by eye. 


r = 0.97. Monte Carlo simulations showed that such high r 
values would only be expected to occur by chance 0.1% of 
the time. 

In the absence of a supernova kick, a positive correlation 
between eccentricity and spin period is expected. Without 
an appreciable kick, the DNS system’s eccentricity follow¬ 
ing the second supernova exp losion is directly prop ortional to 
the mass lost in the process ( Willems et al.|2008l l. A lower- 
mass secondary will evolve niore slowly, prolonging the mass 
transfer phase and allowing the primary more time to accrete 
material and spin up. The low-mass secondary loses relatively 
little mass during its supernova explosion, so the resulting 
DNS system (if it remains bound) has low eccentricity and a 
primary NS with a short spin period. Conversely, a high-mass 
secondary evolves more quickly. A short mass transfer phase 
before the secondary’s supernova leaves the primary with a 
long spin period; greater mass loss during the supernova re¬ 
sults in a DNS system with high eccentricity. 

Using population synthesis simulations of DNS formation, 

S show that this expected relationship 

i ^ 


Dewi et al. 

persists for small supernova kick velocities. The authors as¬ 
sume standard and double-core formation scenarios and that 
most of the matter accreted by the primary NS comes from 
a helium star companion during the He-star-NS binary mass 
transfer phase. Finally, they assume that the accretion rate is 
limited by the Eddington rate for helium accretion and invoke 
a model to account for accretion-induced spin-up of the pri¬ 
mary NS. By applying a small Maxwellian kick velocity dis¬ 
tribution with dispersion a 2 = 20 km s“' for the secondary su¬ 
pernova explosion, the distribution of simulated DNS systems 


et al.j (|2005|l. The Fspin - e relationship is maintained in the 


simulated population for a < 50 km s“*, regardless of the kick 
velocity imparted by the primary supernova. 

We plo t a representation of the simulated DNS popula- 
tion from Dewi et al. ( 2005 1 ) as a shaded region in Figure 
[Keith et al 


.|(|2009[) used a similar representation to argue that 


J1753-2240 is perhaps more representative of the DNS pop¬ 
ulation than previously discovered DNS pulsars, which all lie 
closer to the more sparsely-populated, low-Pspin edge of the 
shaded region. Although J1 930-1852 appears to be approx¬ 
imately consistent with the Dewi et ah] ( 2005l l simulations, 
its location in a sparsely-populated portion of the proposed 
Fspin-e distribution draws the model’s assumptions and con¬ 
clusions into question. 

In their paper, Dewi et al. note that the Apin-e relationship 
is destroyed (randomized) in a simulated population of DNS 
systems with high supernova k ick velocities (a > 5 0 km s“*). 
There is good evidence (e.g. Wong et al. [2010| l that both 
B1913H-16 and B1534H-12 require high supernova kicks, pos¬ 
sibly 3^150 km s“*, but these DNS systems remain com¬ 
pletely consistent with the simulated distribution of systems 
with low kick velocites shown in Figure]^ Therefore, there is 
currently no know n DNS system that supports the claim from 
Dewi et ar[([2005[) that high supernova kick velocity systems 


do not follow a Fspin-e relationship. However, J1930-1852 
(like J1753-2240) provides further evidence that any relation¬ 
ship between spin period and eccentricity is probably much 
broader than was originally thought. 

4.2. Orbital Period-Eccentricity Relationship 

For short orbital period binary systems, one would expect 
a prolonged or more efficient recycling process than for long 
orbital period systems. Since the amount of recycling that oc¬ 
curs is inversely related to both orbital period and spin period, 
these quantities are expected to track one another (i.e. long 
orbital period DNS binaries should have long spin periods). 
We see this relationship in the known DNS population and 
Table 1^ illustrates it nicely. Building on arguments in §4. 1[ 
we expect DNS systems with longer orbital periods to also 
h ave larger eccentricitie s. 

Andrews et al. (2015[l simulate DNS system formation, con- 


sidering three dominant evolutionary channels for the primary 
NS’s companion; (i) wide-orbit common envelope, followed 
by iron core-collapse supernova, (ii) tight-orbit common en¬ 
velope, followed by a second round of mass transfer, then 
iron core-collapse supernova, and (iii) similar to channel ii, 
but a lower-mass secondary He core forms a NS through elec¬ 
tron capture supernova. Using currently understood evolu¬ 
tionary histories for PSRs J0737-3039A/B, B1534H-12 and 
B1913H-16, the authors constrain DNS population models and 
binary parameters. For each synthetic population represent¬ 
ing a given model, they compare confidence intervals that 
result from simulations to the actual DNS population in the 
Pb-e plane. Confidence intervals for all evolutionary chan¬ 
nels combined follow a roughly positive slope in the Pb-e 
plane, as one would expect. High-confidence regions shift de¬ 
pending on the chosen set of model parameters and evolution¬ 
ary channel. For example, higher supernova kick velocities 
tend to produce systems with higher eccentricities; decreasing 
the efficiency of the common envelope ejection mechanism 
produces systems with wider orbits. Of the three evolutionary 
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channels, the first typically has the broadest distribution in or¬ 
bital period. In most cases (i.e. for most chosen sets of model 
parameters), the known DNS population falls well within 3cr 
confidence intervals considering all evolutionary channels (i, 
ii and iii) combined. Comparing J1930-1852’s paramete rs to 
the confidence inte rvals shown for the reference model in I An- 1 
drews et al.|P0151l (Figure 1 in that paper), the combination of 
orbital period ana eccentricity is inconsistent with simulated 
parameters of systems from the three domina nt evolutionary 
channe ls at the 4-5a level. Many models in [Andrews et al.| 
(|2015|l do not produce systems like J1930-1852, which sug- 
gests that it may be useful for constraining theory. 


4.3. Future Work 

PSR J1930-1852’s parameters set it apart from previously- 
studied primary DNS pulsars: the widest orbit (longest orbital 
period), the longest spin period, the largest spin-down rate. 
We hope to determine the pulsar and companion masses inde¬ 
pendently and accurately measure proper motion to infer (i) 
the progenitor mass of the second-born NS and (ii) the mag¬ 
nitude of the s upernova kick it rec eived at birth. Then, fol¬ 
lowing work by jWong et al.| ( |20T0) l, we can make meaningful 
predictions about J1930-1852’s true evolutionary history. 

So far, we detect no Shapiro delay signature in our timing 
residuals (see Figure |^, despite having good orbital cover¬ 
age and a 3.5 hr superior conjunction observation. This non¬ 
detection could be a result of a combination of factors includ¬ 
ing: inadequate timing precision, low companion mass and/or 
small inclination angle. If the non-detection is related to tim¬ 
ing precision, continued timing observations may result in a 


measurement of the delay, and hence an additional constraint 
on the companion mass. 

Upcoming VLBA observations will provide high-precision 
position, proper motion and parallax measurements, which 
will help isolate spin, orbital and post-Keplerian parameters 
in our timing residual fits. Also, because J1930-1852 is rela¬ 
tively nearby (t/oM ~ 1 -5 kpc), we will be able to resolve the 
orbit with VLBA imaging; at the estimated distance, the ma¬ 
jor axis of J1930-1852’s orbit spans ~250 /ias and the VLBA 
routinely obtains 50 /xas precision or b etter with an in-beam 
calibrator available ( jPeller et al.||2012} . Resolving the orbit 
will allow us to place constraints on the system’s inclination 
angle and in turn, the masses of J1930-1852 and its com¬ 
panion. These measurements are critical for eventually un¬ 
derstanding the formation scenario for J1930-1852 and other 
wide-orbit DNS systems. 
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